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Abstract 

There  are  significant  optical  wavefront  perturbationa!  effects  in  the 
use  of  lasers  and  in  imaging  through  atmospheric  turbulence.  In  adaptive 
optics,  an  active  element,  such  as  a deformable  mirror  is  used  to  eliminate 
the  phase  aberrations  of  an  optica:  wavefront.  Diagnostics  to  determine  the 
phase  distortions  are  obtained  by  various  methods.  The  deformable  • bror  is 
used  in  a feedback  loop  in  orde»  to  cancel  the  unwanted  pha  list  >rt ions. 

This  paper  considers  what  is  known  as  a modal  approach  to  fhe  adaptive 
optics  problem  in  which  the  wavefront  distortions  are  expanded  into  1 spatial 
se+  of  polynomials.  The  coefficients  of  the  polynomials  are  temporally 
varying  and  represent  the  time  varying  phenomena logica"1  perturbations  such 
as  atmospheric  turbulence,  thermal  blooming,  and  mirror  distortions.  The 
measurement  devices  for  the  wavefront  are  either  a shearing  interferometer 
or  a Hartmann  array.  The  shearing  interferometer  displaces  or  shears  the 
wavefront  and  then  interferes  the  sheared  imagp  with  the  unsheared  image. 

The  Hartmann  array  yields  a linear  measurement  of  the  aberration  coefficients 
through  an  effective  measurement  of  the  gradient  of  the  wavefron,..  Since  :he 
wavefront  is  temporally  varying,  state  space  models  for  the  aberration 
coefficients  are  obtained.  An  adaptive  estimator  is  developed  in  order  to 
adapt  upon  the  atmospheric  turbulence  structure  constant  and  the  bandwidth 
of  the  atmospheric  turbulence,  and  to  then  obtain  a minimum  mean  square 
estimator  for  the  aberration  coefficients. 

I . Introduction 

In  order  to  obtain  the  maximum  irradiance  of  a laser  beam  on  a distant 
object,  it  is  necessary  to  precisely  control  the  position  of  the  beam  on  the 
object.  The  pointing  control  is  a precision  closed  loop  feedback  system 
which  uses  feedback  measurements  relative  to  where  the  beam  is  currently 
pointing  in  order  to  correct  for  deviations  from  where  it  should  bo  pointed. 
This  is  certainly  an  important  element  in  order  to  obtain  maximum  irradiance 
over  a period  of  time.  However,  without  further  consideration  the  beam 
prop. nation  to  the  object  is  in  essence  an  open  loop  element  in  such  a 
system.  That  is,  there  is  no  feedback  as  to  the  quality  of  the  beam  as  it 
reaches  the  object.  Thus,  the  beam  is  free  to  propagate  according  to  the 
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current  physical  conditions.  Since  the  propagation  is  affected  by  certain 
physical  phenomenon,  the  beam  quality  at  the  object  will  be  reduced 
yielding  a reduction  in  the  irradiance.  These  physical  perturbational 
effects  include  the  initial  laser  wavefront  phase  error,  aberrations  in  the 
optical  train,  atmospheric  turbulence,  and  atmospheric  thermal  blooming. 

These  perturbational  effects  cause  amplitude  losses  and  phase  distortions. 

The  phase  changes  across  the  spatial  extent  of  the  beam  causes  destructive 
interference  at  the  range  of  the  object.  The  irradiance  losses  due  to  the 
phase  distortions  can  be  quite  severe.  Consequently,  it  is  desirable  to 
obtain  feedback  information  as  to  the  beam  quality  at  the  object  such  that 
the  beam  wavefront  exiting  at  the  aperture  can  be  phase  controlled.  The 
phase  control  is  introduced  such  that  at  the  object  range  the  beam  will 
constructively  interfere  to  create  maximum  irradiance.  In  this  way,  the 
beam  is  no  longer  an  open  loop  element  in  the  system  but  becomes  a closed 
loop  subsystem  yielding  all  the  desirable  features  of  closed  loop  control 
such  as  error  correction  and  insensitivity. 

The  technology  area  of  phase  control  of  modification  of  the  wavefront 
exiting  the  aperture  is  known  as  adaptive  optics.  Basically,  this  technology 
requires  a wavefront  diagnostic  subsystem  in  order  to  determine  the  beam 
quality  at  the  object.  The  beam  quality  is  diaqnosed  in  order  to  obtain  a 
measure  of  the  deviation  of  the  wavefront  from  a desired  condition.  The 
controller  subsystem  uses  the  wavefront  diaanostics  to  calculate  the  required 
movement  of  the  controlled  element  in  order  to  obtain  the  necessary  phase 
modifications.  The  controlled  element  can  be  a continuous,  deformable 
mirror  or  segmented  mirrors  that  are  positioned  in  order  to  create  the 
necessary  phase  changes  in  the  outgoing  wavefront.  The  objective,  again, 
is  to  create  the  necessary  phase  changes  in  order  to  obtain  constructive 
interference  at  the  object  range  for  maximum  intensity. 

The  use  of  the  philosophy  of  adaptive  optics  is  not  limited  to  laser 
applications.  It  can  be  used  to  correct  imaginq  system  optics.  This 
correction  is  to  compensate  for  deviations  in  the  object  wavefront  due  to 
perturbational  effects  such  as  atmospheric  turbulence.  Thereby  high 
quality  imaging  may  be  accomplished. 

There  has  been  considerable  effort  in  the  use  of  this  technology  for 
both  laser  and  imaging  applications.  In  references  [ 1 , 2 , 3 [ a segmented 
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mirror  dither  approach  is  considered.  In  particular,  trie  beam  is  split  into 
several  segments  which  are  mdoulated  at  different  frequencies  by  dither 
mirrors.  A detector  measures  the  reflected  intensity  level  modulated  at 
the  different  dither  frequencies.  Synchronous  detectors  were  used  to 
detect  and  obtain  a correction  signal  to  be  added  to  the  mirror  tilt  such 
that  the  field  phasors  become  nearly  aligned  at  the  object  This  condition 
will  yield  maximum  irradiance.  Reference  [4]  considers  the  use  of  estima- 
tion and  control  in  multiditner  adaptive  optics.  In  reference  [5j  a casse- 
grain  telescope  is  autofocused  by  use  of  a sinusoidal  perturbation  t>pe 
adaptive  control  system.  The  focal  plane  f o the  converging  beam  is  adjusted 
in  order  to  coincide  with  the  object  range.  This  yields  [6],  the  maximum 
beam  irradiance  on  the  object.  Reference  [13]  is  devoted  to  the  optical 
technology  of  adaptive  optics.  This  issue  contains  numerous  papers  on 
phase  conjugate  transmitter  systems,  image  sharpening  systems,  transmitter 
muitidither  systems,  and  compensated  imaging.  Another  approach  to  adaptive 
optics  is  the  utilization  of  Bragg  diffraction  [14,15,16]  in  orde1'  to  obtain 
a change  in  the  index  of  refraction  by  d ing  air  with  a sound  wave. 

There  has  been  considerable  interest  and  effort  i u.ing  a closed  loop 
concept  in  optical  systems  for  astronomical  telescopes  [7-12].  Ii  this 
application  the  telescope  mirror  surfaces  ire  disturbed  due  to  structural 
and  thermal  deformations.  A figure  sensor  is  used  to  obtain  the  actual, 
distorted,  mirror  figure.  This  measurement  is  compared  against  the  desired 
figure  and  control  actuators  and  commanded  to  position  the  mirror  surface 
to  minimize  the  error.  In  reference  [11],  a measure  of  image  sharpness  is 
used  as  a performance  index  from  which  a self-optimizing  control  system  is 
used  to  phase  control  the  incoming  wavefront  to  eliminate  phase  distortions. 
See  reference  [13]  for  other  papers  on  image  sharpening. 

The  first  step  in  obtaining  maximum  irradiance  on  an  object  or  for 
image  compensation  is  that  of  wavefront  diagnostics.  This  yields  a meac,:re 
of  the  wavefront  distortion  in  order  to  find  the  correction  signal  to  the 
control  elements.  The  control  elements  can  be  the  adjustment  of  the  focal 
length  of  a r.assegrain  telescope  by  changing  the  distance  between  the 
primary  and  secondary  mirrors  [4,5],  the  control  of  a tilt  mirror,  and/or 
the  rontro1  of  a deformable  m’rror. 
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The  measurements  of  the  wavefront  distortions  are  noisy.  In  particular, 
detector  noise,  background  noise,  and  photon  noise  for  photon  limited  appli- 
cations limit  the  accuracy  of  measurement  of  the  distortions.  This  paper 
develops  the  optimal  estimator  for  estimating  the  wavefront  distortions  with 
two  types  of  measurement  devices,  a shearing  interferometer  and  a Hartmann 
array.  Although  these  are  two  physically  different  optical  devices,  the 
mathematics  of  the  measurements  are  similar.  The  estimators  are,  thus, 
similar  and  it  is  reasonable  to  consider  both  in  the  same  paper. 

The  distortions  are  expanded  into  an  orthonormal  set  of  polynomials 
where  the  coefficients  are  called  aberration  coefficients.  This  corresponds 
to  the  proper  optics  terminology  and  the  coefficients  correspond  to  well 
known  optical  errors  such  as  focus,  astigmatism,  coma,  etc.  The  aberration 
coefficient  temporal  properties  are  modeled  as  linear  Markov  models.  The 
bandwidths  of  the  models  are  directly  proportional  to  the  relative  wind 
speed  across  the  optical  aperture.  The  driving  noise  variance  is  directly 
proportional  to  the  atmospheric  turbulence  structure  constant,  Cn:'  (see 
reference  [17]). 

This  paper  is  the  first  application  of  adaptive  estimation  to  adaptive 
optics.  The  adaptive  estimator  converges  readily  to  the  required  estimates 
of  the  aberration  coefficients  in  the  presence  of  uncertain  atmospheric 
turbulence.  Thus,  the  adaptive  estimation  of  aberration  coefficients  in 
adaptive  optics  is  an  extremely  useful  and  important  technique  to  enhance 
performance  of  high  energy  lasers  and  to  enhance  performance  of  imaging 
systems.  The  adaptive  optics  application  represents  a statn  of  the  art 
area  in  optics.  Thus,  the  paper  contributes  also  to  the  optical  literature. 

The  paper  is  broken  into  six  sections.  Section  II  contains  the  problem 
statement.  Section  III  gives  the  measurement  devices  and  their  equations. 
Section  IV  considers  the  dynamic  models  for  the  aberrations.  Section  V 
contains  the  formal ation  of  the  adaptive  estimator,  and  Section  VI  contains 
the  application  of  the  estimator  to  the  shearing  interferometer  problem. 

1 1 . Problem  Statement 

In  order  to  correct  the  wavefront  of  an  optical  beam,  it  is  necessary 
to  spatially  modify  the  phase  of  the  beam.  There  are  several  ways  of 
accomplishing  a phase  change,  i.e.,  phase  plates  f 1 8 1 , Bragg  diffraction. 
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focal  length  and  tilt  adjustment,  and/or  deformable  mirrors.  There  is  a 

need  to  change  the  phase  in  real-time  in  order  to  compensate  for  temporally 

as  well  as  spatially  changing  distortions,  and  do  this  with  low  power 

requirements.  In  high  energy  laser  applications,  the  most  viable  methods 

are  that  of  focal  length  ancj  tilt  adjustment  with  higher  order  aberration 

/ 

correction  using  a deformable  mirror.  The  mirrored  surfaces  must  be  cooled, 
however. 

The  intensity  of  a laser  beam  at  the  far  field  is  a mapping  of  the 
phase  as  well  as  the  amplitude  [18,19|.  A ratio  of  the  ideal  intensity 
without  phase  distortions,  I*,  to  the  intensity  with  phase  distortions,  I, 
is  called  the  Strehl  ratio,  i.e.. 


It  may  be  shown  | 1 8 . that  the  Strehl  ratio  is  always  less  than  one  in  the 
presence  of  phase  distortions.  High  quality  imagery  also  requires  that  the 
jhase  ' t e z The  • ->c.  : • • • /mb  • 1 i tally  in  figure  1 . 

In  both  the  problems  of  imaging  and  high  energy  lasers,  atmospheric 
turbulence  acts  to  phase  distort  the  optics.  In  addition,  base  vibrations 
cause  effective  phase  distortions,  thermal  blooming  of  high  energy  lasers 
are  a problem,  and  even  the  resonance  in  the  laser  cavity  cause  distortion 
problems.  However,  if  the  phase  distortions,  $(x,y),were  known  over  all 
the  spatial  regime,  then  a deformable  mirror  along  with  focus  and  tilt 
adjustment  can  be  used  to  eliminate  much  of  the  distortion.  In  order  to 
accomplish  + his  diagnostic  problem,  the  wavefront  distortions,  I(x,y)  may 
be  expanded  into  a set  of  spatial  polynomials  defined  over  the  spatial 
extent  of  the  aperture.  In  particular,  the  polynomials  are  chosen  to  be 
orthogonal  over  the  aperture.  Two  sets  are  common  in  optics  work.  Thp 
first  set  of  polynomial  expansions  is  known  as  Zernike  polynomials  118). 

The  second  set  is  used  originally  by  Fried  |17j  and  later  by  Hogge  [?0|  in 
atmospheric  turbulence  work.  The  second  expansion  can  be  defined  as  follows 


•' ( x ,y ) = £ a - F -(x,y) 


(2) 
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where  the  domain  of  definition  for  x and  y is  the  circular  region  of  the 
aperture.  The  a-'s  are  the  aberration  coefficients  and  F.  are  a set  of 
orthonormal  polynomials  related  to  the  Zernike  polynomial.  The  first  ten 
polynomials  and  their  physical  significance  are  given  below. 


Fj  (x,y) 
F?(x,y) 
F,(x,y) 


1 


(uniform  phase  shift  across  the  aperture) 


(tilt  across  aperture) 


Fu(*.y) 

F5(x,y) 

F.(x,y) 

6 

Fy(x,y) 

Fs(x,y) 


Fg(x,y) 


Flc(x,y) 


(refocus) 


(astigmatism) 


(coma) 


(3) 
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are  the  dominant  terms.  Furthermore,  due  to  the  fact  fhat  other  physical 
phenomenon  such  as  thermal  blooming  and  mirror  thermal  distortions  are  also 
adequately  represented  by  the  low  order  aberrations,  the  wavefront  distor- 
tions may  be  adequately  assumed  to  be  represented  as  a finite  polynomial 


*(x,y)  = 5T  a.F.(x,y).  (4) 

j=l  J J 

In  atmospheric  turbulence,  the  aberration  coefficients  are  not  constant 
Rather,  they  are  temporally  varying  with  bandwidth  related  to  the  relative 
wind  speed  across  the  aperture  and  with  variance  that  is  a function  of  the 
turbulence  structure  constant,  C- , which  is  a measure  of  the  turbulence 
strength.  Thus,  it  is  more  approoriate  to  write  equation  (4)  as 


n 

( x , y , t ) - X/  a-(t)F  (x,y).  (5) 

j=l  J J 

Given  the  ph.ve-  distortions  as  in  equation  (5),  the  deformable  mirror  may  be 
controlled  to  have  a phase  eaual  and  opposite  the  phase  distortions  such 
that  the  total  phase  is  zero.  This  leads  to  another  reason,  less  physical 
than  the  dominance  of  the  lower  order  terms,  for  the  expansion  of  c into  a 
finite  po  . o<  ' ■ l . That  is,  the  control  of  the  deformable  mirror  is  limited 
to  the  highest  spatial  frequencies  (as  well  as  temporal  bandwidth)  that  may 
be  controller  due  to  the  finite  spacing  o control  actuators  on  the  mirror. 

In  atmospheric  turbulence,  the  image  of  the  far  field  may  b^  used  to 
determine  the  chase  distortions  from  which  the  reciprocity  principle  may  be 
used  to  then  a-  el  the  phase  via  the  deformable  mirror  control.  Systems 
which  use  th  . t’.pe  of  wavefront  measurement  are  known  as  imaqe  compensating 
systems  |li].  A measurement  of  the  phase  distortions  from  the  image  may  be 
obtained  by  using  a shearing  interferometer  or  by  a Hartmann  array.  The 
measurement-,  of  these  device  may  be  related  to  the  aberration  coeff indents 
plus  noise.  The  dynamic  model  (or  the  aberrations  are  uncertain  models  in 
that  they  depend  upon  knowledge  of  the  relative  wind  speed  as  well  as  the 
turbulence  structure  constant.  Since  the  measurements  are  noisy  and  the 
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dynamic  models  contain  uncertain  parameters,  it  is  necessary  to  use  adaptive 
estimation  in  order  to  obtain  an  estimator  for  the  coefficients.  An  impor- 
tant aspect  of  applications  of  the  estimator  is  that  the  convergence  times 
required  are  extremely  fast.  Thus,  the  adaptive  estimation  scheme  chosen 
includes  a probability  estimator  for  the  uncertain  coefficients  although  the 
computational  burden  is  greater  than  other  adaptive  estimation  schemes,  the 
convergence  times  required  justify  the  use  of  this  algorithm.  The  next 
section  gives  a basic  description  of  the  measurement  devices  as  well  as 
their  equations. 


III.  Measurement  Devices 

There  are  two  major  measurement  devices  used  in  image  compensation 
systems.  The  first  is  that  of  a shearing  interferometer  and  the  second  is 
that  of  a Hartmann  array.  The  measurement  of  the  wavefront  using  a shearing 
interferometer  is  depicted  in  Figure  2 for  a one  directional  shear.  The 
wavefront  is  sheared  in  two  orthogonal  directions.  The  interference  wave- 
front  in  one  direction  may  be  written  as 

V(x,y;sx,t)  = :(x  + sx,y,t)  - t(x,y,t)  (6) 


where  j>  is  the  phase  front  under  test,  sx  is  the  lateral  shear  direction, 
and  V is  the  quantity  to  be  measured. 

Detector  arrays  behind  the  wavefront  interferometer  using  zero  crossing 
will  measure  the  wavefront  plus  noise.  Thus  (assuming  that  two  detector 
arrays  are  available  to  measure  the  vertical  and  the  horizontal  shear),  the 
measurements  for  the  i-th  detector  in  each  array  is  given  as 


V.  . ( x ,y ; 


Sj.t)  . 


. .(x+s-iS  .,y+s-''  -,t)  - <!> ( x . ,t) 
ijv  J xj  J .1  yj  iy 


i ] 


(7) 


i - 1,2 q 

j = x,y 


where  . and  - . are  zero  when  x f j or  y t j respectively,  and  when  i 

* j .yj 

represents  the  i-th  detector  in  the  j-th  array  (q  detectors  in  each  array), 
j = x corresponds  to  a horizontal  shear,  j = y corresponds  to  a vertical 


shear,  and 


ij 


is  the  detector  noise  on  the  i -j-th  detector.  The  noise 


ij 
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is  assumed  to  be  zero  mean,  white  with  covariance  r-  . it  is  assumed  that 

i J 

the  detector  noise  terms  are  independent  of  other  detectors.  The  shearing 
distance  is  chosen  equal  to  the  separation  of  the  detectors  i.>  the  array. 
The  wavefront  measurements  may  be  placed  into  an  array 

y = Tt  + n (8) 


where  y is  a vector 


T 

y 


2 ,x’ 


q,x 


2,y 


is  a vector  formed  from  the  independent  phases  in  equation  (7),  and  n is 
a noise  vector,  zero  mean,  white  with  diagonal  covariance 


E|  (t),  (0Tj  = Ro  ( t-  r ) 

where  R is  diagonal  with  elements  r-^.  The  measurements  in  equation  (7) 
will  be  such  that  some  of  the  and  . will  correspond  to  the  same  spatial 

phases . Thus , there  will  only  be  a fewer  number  of  phases  than  4q  phases. 

Th  matrix  T -' $ a matrix  of  plus  and  minus  ones  and  zeros.  A typical 
example  is  shown  in  Figure  3 in  which  the  shear  measurements  may  be  written 
in  terms  of  the  phase  values  at  the  i — th  spatial  position.  The  12  measure- 
ments n vector  form  a-’e 
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The  wavefront  distortions  may  be  expanded  into  tne  polynomials  as  in  equa- 
tion (5) , i .e  . , 


(V.yrt)  = Ij  a.(t)Fj(x.,yi) 


(10) 


where  (x.^. ) corresponds  to  the  i-th  spatial  point.  Thus,  a vector  equa- 
tion may  be  written  for  the  phases  in  terms  of  the  aberration  coefficients, 
i .e.  , 


F (Xj ,y, ) , Fr(Xj ,y} ),  ....  Fn(x,  ,y, ) 

(t) 

F;(x?,y?),  F (x  ,yj Fn(x  ,y  ) 

a (t) 

A 

3 

Fj  (x3,y3) , F.,  (x3  ,y3 ) , ....  Fn(x3,yJ 

a , ( t ) 

A 

' m 

F i ^xm,ym^  ’ r (xm’ym) 'VvV 

am(l)_ 

(ID 


or  equation  (11)  may  be  written  in  vector  form  as 
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__ 


4'  = F a ( t ) . 


02) 


Equation  (12)  may  be  used  with  equation  (8)  to  obtain  a measurement  equa- 
tion in  terms  of  the  aberration  coefficients,  i.e.. 


y = ( TF ) a + o 
= Hsla  + n 


where  the  definition  of  is  obvious  and  the  subscript,  SI,  denotes  the 
measurements  using  a shearing  interferometer.  This  is  the  desired  measure- 
ment  equation  for  the  shearing  interferometer. 

The  measurement  of  the  wavefront  using  a Hartmann  array  is  depicted  in 
Figure  4.  Basically,  the  image  is  formed  on  the  Hartmann  array  by  a lens 
system.  Since  the  distance  to  the  object  is  long,  a wavefront  impinging  on 
the  array  would  be  approximately  planar  without  turbulence  effects.  How- 
ever, turbulence  distorts  the  phase  of  the  wavefront  as  in  Figure  1.  Thus, 
the  image  across  the  Hartmann  array  is  phase  distorted  in  the  presence  of 
atmospheric  turbulence.  The  Hartmann  array  consists  of  several  hundred 
pinholes  (or  fly  eye  lenses)  that  deflect  the  light  back  to  quad  detectors 
located  behind  each  hole  (or  lens).  The  quad  detectors  will  yield  measure- 
ments proportional  to  the  gradient  of  the  wavefront  in  both  x and  y direc- 
tions, i.e.,  for  each  pinhole  two  measurements  are  obtained,  i.e., 


where  q is  the  number  of  pinhole-detector  pairs.  The  gradient  of  the  wave- 
front  in  terms  of  the  aberration  coefficients  may  be  obtained  by  differen- 
tiating equtMon  (5),  i.e.. 
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Equations  (14)  and  (15)  may  be  put  into  vector  form  as 
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or  in  vector  form 

y = H a + n (19) 

where  y is  the  measurement  vector  in  equation  (18),  is  the  measurement 
matrix  of  partials,  and  n is  the  vector  of  measurement  noises.  The  measure- 
ment noise  is  assumed  zero-mean,  white  with  covariance 


E |i( t) n( t ) j = R(t)S(t-i) 


where  R is  in  block  diagonal  form  since  the  measurement  noise  nx.  is 

correlated  in  general  with  nv  , v-  = l,2,...,q.  This  particular  aspect 
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complicates  the  estimator  somewhat  as  each  pair  of  measurements  for  one 
detector  must  be  processed  through  the  estimator  rather  than  sequentially 
processing  each  measurement.  However,  the  form  of  the  estimator  is  still 
simplified  over  the  form  required  if  R was  not  block  diagonal. 

This  section  then  contains  the  development  of  the  measurement  equa- 
tions for  both  the  shearing  interferometer  and  the  Hartmann  array  image 
compensation  systems.  The  equations  are  in  the  form  of  linear  transforma- 
tions of  the  aberration  coefficients.  This  allows  for  the  development  of 
the  instimation  equations  for  both  systems.  The  next  section  considers 
approx imations  for  the  state  space  models  of  the  aberration  coefficients. 

I V . Dynam i c Model  ing  of  jth_e  Aberration  Coefficients 

In  order  to  develop  an  estimator  structure""  for  the  estimation  of  the 
aberration  coefficients,  it  is  necessary  to  obtain  models  for  the  temporal 
variation  of  the  coefficients.  The  models  must  be  such  tnat  they  adequately 
represent  the  actua7  statistics  of  variation,  and  yet  be  simple  enough  for 
implementation  into  the  filter.  Furthermore,  the  models  will  be  shown  to 
be  functions  of  the  atmospheric  turbulence  bandwidth  as  well  as  the  atmos- 
pheric turbulence  structure  constant,  Cn2  (see  references  [24,  25]  for  a 
precise  definition).  These  models  are  approximations  similar  to  that  of 
reference  f 26 1 for  the  intensity  fluctuations  of  a laser  beam  arid  similar 
to  reference  [27]  for  the  temporal  phase  distortions.  However,  since  the 
estimator  is  adaptive,  it  will  adapt  upon  the  best  model  structure  for  the 
particular  realization  of  atmospheric  turbulence.  Also,  the  time  interval 
for  use  of  a high  power  laser  as  well  as  imaging  through  atmospheric  turbu- 
lence is  short.  Thus,  the  models  are  adequate  for  this  problem  and  will 
yield  an  estimator  structure  that  can  be  implemented  real  time. 

In  Fried's  paper  fl7J  on  wavefront  distortion,  it  is  shown  that  the 
mean  square  value  of  the  linear  tilt  coefficients  is  given  as 

E(aL7)  = 0.694D2(D/rQ)5/3  (20) 

where  a^-  is  equal  to  (a,'  +a  ),  the  mean  square  value  of  the  spherical 
coefficients  is  given  as 

5/3 


E(as2)  = 0.01657D2(D/ro) 


(21) 


where  au-  is  equal  to  a^- , and  the  mean  square  value  of  the  quadratic  coeffi- 
cients is  given  as 

EU^7)  = O.Ob27D(D/ro)5/3  (22) 

where  aq  is  equal  to  (a.'  t-ar:  +a(:).  The  diameter  of  the  aperture  is 
denoted  as  D,  and  the  quantity  rp  is  a length  parameter  defined  as 

rQ  = (6.88/a)3/5  (23) 

where 

a = 2 91  ( 2 •:  / > ) 2 f C 2(s)Q(s)ds  (24) 

"'path 


where  \ is  the  wavelength  and  where  Q(D)  is  a weighting  factor  depending  on 
the  nature  of  the  optical  source.  In  the  case  of  an  infinite  plane  wave 
source  with  a constant  structure  constant,  the  parameter  rQ  is  given  as 


6.88 

2.91  k' RC  • , 
n -J 


3/5 


meter s 


(25) 


where  k = 2n/A  is  the  wavenumber,  R is  the  range  in  meters,  and  Cn2  is  the 
turbulence  structure  constant  which  is  a measure  of  turbulence  strength. 
References  [20,  28]  discuss  the  temporal  variation  of  the  atmospheric  turbu- 
lence. It  is  shown  that  the  corner  frequency  of  the  turbulence  for  each 
aberration  is  of  the  order  of  the  relative  wind  speed  divided  by  the  diameter 
of  the  aperture.  In  reference  f 28 ) approximations  to  the  power  spectral 
dens' ties  are  given.  It  is  assumed  that  adequate  approximations  to  these 
models  for  each  aberration  are  of  the  form 

a • = - 6 ^ a i + v'YpT.  o^,  (26) 

i-2,3, . . . ,n. 


where  the  ■ is  given  as 


(27) 
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with  k.j  given  in  Table  I,  where  V is  the  relative  wind  velocity,  and  where 
a.j  is  zero  mean  white  noise  with  unity  variance.  The  variance  may  be 
calculated  by  use  of  Fried's  paper.  Thus,  the  terms  are 


and 


i=2,3. 


°4  = ^(V), 


i=5,6. 


(28) 


It  is  assumed  in  such  a set  of  models  that  the  aberrations  are  uncorrelated. 
They  are  not  in  the  expansion  functions  as  given  in  [20].  However, 
uncorrelated  expansion  functions  may  be  found.  Also,  additional  uncertain 
model  parameters  may  be  used  to  ascertain  this  correlation.  These  parameters 
may  be  augmented  to  the  adaotive  estimator. 

The  first  coefficient  which  appears  as  a constant  phase  shift  doesn't 
affect  the  far  field  intensity.  It  may  be  noted  that  the  variances,  o., 
are  functions  of  the  turbulence  structure  constant.  Since  the  turbulence 
strength  is  not  known  a priori  for  a ptach  of  turbulence,  the  structure 
constant,  Cp2,  is  uncertain.  Furthermore,  since  the  relative  wind  speed  is 
not  known  precisely  for  an  aircraft  scenario,  V may  be  uncertain.  Thus, 
the  estimator  must  adapt  on  the  uncertain  parameters  in  order  to  learn  them 
for  use  in  the  estimation  processes.  Furthermore,  structure  adaptation  may 
be  used  to  refine  the  models  for  a particular  realization  if  the  time  inter- 
val demands  a more  accurate  model . 

The  model  approximations  are  reasonable  for  the  logic  previously  stated. 
That  is,  other  work  in  laser  fluctuations  have  used  similar  models,  the 
estimator  will  adapt  upon  uncertain  parameters,  the  analytics  for  turbulence 
may  not  correspond  to  a given  scenario,  and  a computationally  feasible  model 
is  necessary  for  real-time  implementation. 

The  next  section  gives  the  structure  of  the  adaptive  estimator. 
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V • Adaptive  Estimator 


This  section  contains  the  equations  for  the  adaptive  estimator.  The 
estimator  is  based  upon  the  original  work  by  Magi T 1 [29|  and  extended  by 
Lainiotis  |30|  and  Hawkes  and  Moore  [31].  In  particular,  as  the  measure- 
ments for  the  shearing  interferometer  are  in  a form  whereby  the  measurement 
noise  covariance  matrix  is  diagonal  and  as  the  measurements  for  the  Hartmans 
array  are  in  a form  whereby  the  measurement  noise  covariance  matrix  is 
block  diagonal  with  blocks  of  2x2,  the  sequential  processing  adaptive 
estimation  form  of  Hawkes  and  Moore  is  especially  viable.  In  this  form, 
the  measurements  may  be  processed  in  a sequential  form  (one  at  a time).  This 
only  requires  a scalar  inverse  for  the  shearing  interferometer  measurement 
and  a 2x2  inverse  for  the  Hartmann  array. 

The  dynamic  equation  for  the  aberration,  equation  (26),  may  be  dis- 
cretized as 


a.(k+l)  = exp  ( - 1-: i T ) a i ( k ) + Qi(k)u.(k)  (29) 

whe  re  , 

/t(a.') " 

qi(k)='/'2ir_  t1  ■ exp  ("2  BiT)] 

u-  is  zero  mean  white  noise  with  unity  variance,  and  T is  the  sampling 
period.  The  dynamic  equations  may,  thus,  be  written  in  the  general  form 


a(  k+1 ) = ■!>  ( v ) a ( k ) + r(v,Cn2)u(k)  (30) 


where 


J(v)  = 


exp(-k  jj) 
0 
0 
0 
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0 

exp(-k, 

0 

0 

0 


0 

0 

exp(-klt  jj) 
0 
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exp(-k.  ~) 
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0 
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exp(-k  q) 
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r 


and  where 


(v’Cn2)  - 


q.(k) 

c. 

0 

0 

0 

0 


q3(k) 

0 

0 

0 


0 

0 

q,  (k) 

4 

o 

0 


0 

0 

0 


q5(k) 


0 

0 

0 

0 

q,(k) 


where 


q n-  ( k ) = qi(k,v,Cn;) 


from  the  definition  of  q^k).  It  is  assumed  that  v and  Cn  are  defined  over 
a discrete  range.  The  two  parameters  may  be  contained  in  the  finite  set 


ve;vi’v2 vq } [ 


and 


r 2-  lr  2 r 2 
n E jl'n1  ,Ln?  > • • • » 


with  a priori  probability  density  function 

q, 


C 2 * 

V 


P(v)  = E F* r ( v i ) 6 ( v-V i 


i=l 


and 


P(C 


n . pr(cn;n(cn'-cnp. 


The  a priori  density  function  for  the  pair  is  given  (assuming  independence) 
as 

q,  q, 

(31) 

t r • l ■ r • "n  • i " i 


q.  q. 

P(v,Cn2)  = E E Pr(vi)Pr(Cnj2)6(v-vi,Cn2-Cn2) 
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or  by  defining  = {v,Cp- 


P(e)  = £ Pr('  ,- ) • ( -i  •) 

M r 1 J 


(32) 


where  . is  defined  for  each  q,  x q,  pair  of  v’s  and  Cn  ‘s  and  q = q.  xq 
The  algorithm  for  the  shearing  interferometer  measurements  will  be 
written  explicitly  for  completeness.  The  Hartmann  array  algorithm  is  an 
immediate  extension  of  this  and  will,  thus,  not  be  written. 

Equation  (13)  may  be  partitioned  as 


Hc , 


q 


(33) 


hi  r<  - • is  t r m of  the  matrix.  Tne  e ■ ement  equation  may  thus 

be  written  as 


y; 

"hj  -3  + S " 

II 

h?  a + n. 

Uc  • • 

0 

h.  a + r 

L C’J 

q q 

(34) 


Using  the  results  in  [30],  the  optimal  estimator  for  the  aberration 
coefficients  given  the  measurements  up  to  and  including  time  k is 


a+(kU'k)  = iC  a+(k' ik,ei)Pr(eii  ik) 

i = l 


(35) 


where  a4(k  t'k*e.j)  is  the  fK -condi ti onal  estimate  for  a and  Pr(a^ | k)  is  the 
probability  of  e given  the  measurements.  The  covariance  of  this  estimate 
is 
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— -n-  - — ■ 


P 


q 

p(k)  = iC{p(k|  °-j)  + [a(kUk,e.)  - a+(k|tpk)][a+(k|ipk,ei)  - a+(k|\|»k)]Tj . 
The  ©^-conditional  estimates  may  be  found  between  measurements  as 

a"(k+l|*k,ei)  = 4>(0i)a+(k|n>k,e.), 

i = l ,2, . . . ,q 
and 

p-(k+1'e.)  = Mei)p+(k!e.)$T(0.)  + r(0.)r(0.)T, 

i = l ,2,. . . ,q 

where  $>(o.)  and  r(0.)  are  as  in  equation  (30)  evaluated  with  the  i-th  6 
a measurement  the  n ..-conditional  estimate  may  be  found  by  sequentially 
processing  the  measurements  through  the  following  algorithm.  The  first 

measurement  processing  is 

al  i(k  ’kVr(,i)  = a"(k!^k_1  ,ei)  + Kk(l)(01)[y1-h1a-(k|*k.1,01)] 

i=l ,2 , . . . ,2q 

where 

Kk(l)(oi)  = P'(k|  j)h1T[h1P'(k|0i)h]T  + r , ] ” 1 , 

i = l ,2, . . . ,2q 

where  r,  is  the  covariance  of  the  final  measurement  noise  element  and 
a ^ ^ ( k ! ^k_ -j  »y j » ■ ’ - ) denotes  the  estimate  given  the  past  measurements  and 
measurement  on  the  first  row  of  the  measurement  vector.  The  covariance 
this  estimate  is 


p(!)(k|o.)  = [l-Kk{l)(ei)hI]P'(k!ei), 

i=l ,2, . . . ,2q 

The  required  estimate  is  given  by  cycling  through  the  iterations 


(36) 

(37) 

(38) 

. At 

(39) 

(40) 

the 

of 

(41) 
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a('  )(!<:  .k,y.,y  ,, y?.ei)  - a('  '(kUk,y1  ,y2, y?_5 J + 


Kk!  ;1V[A- 


iui":1(i!*jt.y1  ,y7,- ■ • .y  ■«;)] 


7=2,3, .. . ,2q 


where 

K,(,)(  ■)  = P(?'l)(k!ei)h?I[h;P(f'l)(k|ei)  . h.T+  r ]_1  (43) 

where  r is  the  covariance  of  the  f.-th  measurement  noise  element.  The 
covariance  update  is  given  by 

P(,)(k  i)  = [I-Kk('  )(ei)hjp(f'i)(kjei),  (44) 

7=2,3, . . . ,2q 

The  required  --conditional  estimate  is  given  by 

a+  ( k | k , e i ) = a(2q)(k;  vk,  i)  (45) 

and  its  covariance  is  given  as 


P+(k!  .)  = P^2q^(k|ei). 


It  may  be  noted  that  only  a scalar  inverse  need  be  taken  to  find  the  required 

estimate.  In  order  to  update  the  necessary  probabilities  sequentially,  it 

(O  T ~ (z) 

is  necessary  to  store  [hfPv  ;(k|e^)h;  + r. ] and  the  residuals  y . - h - a' ' ' (k j 9 . ) 

for  7=l,2,...,2q  and  i = l,2 q.  These  quantities  are  available  as  they 

are  calculated  from  the  previous  equations.  The  required  probabilities  are 

then  calculated  as  in  reference  [31],  i.e.. 


VV  !k>  * "T1 


(47) 


for  i=l ,2, ... ,q  where 


c(  '^i)  = yz ' h4a(e)(kK|c_1.y1,y2,...,yJl_1.0i) 


(48) 


and 


^k(0(ei)  = Kpl*;(kiej)V  + rJ- 


id, 


(49) 


The  terms  v^^(  • ) and  £2.^(6,.)  have  been  stored  during  the  previous  calcu- 


lations. 


VI . Results  and  Conclusions 

The  estimator  was  simulated  for  various  types  of  lasers,  several  different 
propagation  lengths,  and  several  optical  aperture  diameters.  The  three  types 
of  lasers  simulated  were  the  aas  dynamic  laser  (6DL)  at  10.6  x 10"1  meter 
wavelength,  the  electric  discharge  laser  (EDL)  at  9.28  x 10"6  meter  wave- 
length, and  a chemical  laser  using  deuterium  floride  (DF)  at  3.69  x 10“6 
meter  wavelength.  Both  short  ranges  1-5  kilometers  and  medium  ranges  20-50 
kilometers  were  used.  Several  different  signal  to  noise  ratios  were  used. 

All  units  yield  a phase  error  in  radians.  The  tilt  errors  are,  for 
example,  in  units  of  radians/meters.  Since  the  resolution  of  any  plots 
would  be  less  than  minimal  because  of  the  significant  decrease  in  several 
orders  of  magnitude  of  the  covariance  of  estimation  error,  plots  of  the 
results  will  not  be  shown.  Instead,  both  the  initial  and  steady  state 
covariance  of  the  estimation  error  for  each  case  will  be  given  for  each 
aberration  order.  The  steady  state  solution  is  a strong  function  of  the 
measurement  noise  and  less  of  a function  of  the  dynamic  models.  The  cases 
shown  are  for  an  approximate  signal  to  noise  ratio  of  five.  The  sample  rate 
for  the  measurements  is  0.0005  seconds.  This  is  considered  realistic  with  a 
state  of  the  art  computer.  Five  possible  values  of  Cp2  and  the  bandwidths 
were  used.  The  actual  Cn2  for  each  run  was  1 x 10-16.  The  turbulence  band- 
widths  were  chosen  with  a transonic  relative  wind  speed. 

For  the  case  of  a GDL  with  one  meter  optics  and  a range  of  one  kilometer, 
the  initial  standard  deviation  for  the  two  tilts,  refocus,  and  two  astigma- 
tisms were  0.7179,  1.56  x 10*1,  and  1.338  x 10“*,  respectively.  In  several 
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Monte  Carlo  runs,  the  correct  parameters  were  learned  before  0.02  seconds. 

The  steady  state  standard  deviation  at  0.025  seconds  after  a measurement  were 
as  follows:  x-tilt,  2.43  x 10  ' ; y-tilt,  2.42  x 10"3;  refocus,  7.14  x 10"4; 

quadratic  astigmatism,  6.86  x IQ*4;  and  product  astigmatism,  7.74  x 10"4. 

The  x-tilt,  for  example,  would  correspond  to  an  uncertainty  of  1.52  x 10"-: 
wavelengths.  With  an  initial  range  of  50  kilometers,  the  initial  standard 
deviations  for  the  two  tilts,  refocus,  and  two  astigmatisms  were  5.259,  1.11, 
and  0.946,  respectively.  Again,  the  true  parameters  were  learned  before  0.02 
seconds  in  several  runs.  The  steady  state  standard  deviations  at  0.025 
seconds  afte^  a measurement  were  as  follows:  x-tilt,  4.665  x 10'  y-tilt, 

4.47  x 10" refocus,  2.93  x IQ'1;  quadratic  astigmatism,  3.22  x 10'  , and 
product  astigmatism,  2.687  x 10" It  can  be  seen  that  the  increase  in 
propagation  distance  increases  the  initial  and  the  final  uncertainty  The 
tilt  uncertainties  at  steady  state  are  about  doubled  for  an  increase  of  range 
by  50  ti^tes. 

For  the  case  of  a GDL  with  0.6  meter  optics  and  a range  of  5 kilometers, 
the  initial  standard  deviations  for  the  two  tilts,  refocus,  and  two  astigma- 
tisms were  1.605,  3.51  x 10“’,  and  2.q9  x 10"',  respectively.  In  several 
Monte  Carlo  runs,  the  correct  para  vere  irm  * before  0.02  seconds. 

The  steady  state  standard  deviations  at  0.9.5  seconds  after  a measurement 
were  as  fellows:  x-tilt,  3.46  x 10”  , y-tilt,  3.38  x 1 0'  1 ; refocus,  1.415  x 

10"';  auadratic  astigmatism,  1.419  x 10"  and  product  astigmatism,  1.413  x 
10"  \ 

For  the  case  of  an  EDI  with  0.6  meter  optics  and  a range  of  5 kilometers, 
the  initial  standard  deviations  for  the  two  tilts,  refocus,  and  two  astigma- 
tisms were  0.'187,  1.57  x 10“’-,  and  1.339  x 10'1,  resDecti vely . The  correct 
parameters  were  lean  d b°fore  0.035  seconds.  The  steady  state  standard 
deviations  at  0.950  seconds  after  a measurement  were  as  follows:  x-tilt, 

2.43  x 10"  ; y-tilt,  2.426  x 10"’;  refocus,  7.15  x 1 0" 4 ; quadratic  astigma- 
tism, 6.87  x 10"';  and  product  astigmatism,  7.74  x 10"4. 

For  the  case  of  a OF  laser  with  0.6  meter  optics  and  a range  of  5 kilo- 
meters, the  in  Hal  standard  deviations  for  the  two  tilts,  refocus,  and  two 
astigmatisms  were  4.61,  1.008,  and  8.59  x 10";,  respectively.  The  steady 
state  standard  deviations  at  0.02  seconds  after  a measurement  were  as 
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follows:  x-tilt,  4.60  x 10*3;  y-tilt,  4.^1  x 10*  ; refocus,  2.79  x 10  ' , 

quadratic  astigmatism,  3.045  x 10*3;  product  astigmatism,  2.569  x 10*'. 

\ * 

With  a range  of  50  kilometers  and  1 meter  optics  the  initial  standard  devia- 
tions are  respectively  14.58,  3.19,  and  2.72.  At  steady  sta*e,  the  standard 
deviations  are  as  follows:  x-tilt,  5.11  x 10*3;  y-tilt.  4.85  x iO*  , 

refocus,  4.39  x 10*  ';  quadratic  astigmatism,  5.338  x 10* and  product 
astigmatism,  3.93  x 10"3. 

It  may  be  seen  that,  as  is  well  known,  the  uncertainties  for  shot  er  wave- 
lengths through  atmospheric  turbulence  are  initially  greatly  increased. 

However,  the  adaptive  estimator  is  able  to  decrease  this  uncertainty,  which 
is  much  higher  than  the  uncertainties  at  larger  wavelengths,  to  a -atue 
nearly  the  same  as  the  larger  wavelength  lasers.  The  estimator  is  seen  T.o 
bring  the  uncertainty  levels  down  to  that  required  for  control  while  adapting 
upon  the  uncertain  atmospheric  turbulence  strength  and  bandwidth  Other  runs 
were  accomplished  at  lower  and  higher  signal  to  noise  levels  which  indicate 
similar  results. 

This  paper  introduces  the  adaptive  optics  problem  to  the  control  community. 
The  adaptive  estimation  structure  developed  within  is  a viable  structure  for 
estimation  of  aberration  coefficients  in  adaptive  optics  systems.  The 
estimator  has  been  well  exercised  with  several  Monte  Carlo  runs  for  GDL, 

EDL,  and  DF  lasers  at  several  ranges  and  several  different  optical  apertures. 
The  estimator  may  be  implemented  into  a minicomputer  with  ether  time  varying 
or  constant  gains  for  real  time  applications. 


UNCORRECTED  LASER  BEAM  PROPAGATING  THROUGH  ATMOSPHERIC  TURBULENCE 
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FIGURE  1 LASER  AND  IMAGE  PROPAGATION  THROUGH  ATMOSPHERIC  TURBULENCE 


PHASE  NATURE 


FIGURE  3 TWO  DIMENSIONAL  SHEARING  INTERFEROMETER 
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